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ABSTRACT 

It is suspected that many strong gravitational lens galaxies lie in poor groups or rich clusters of galaxies, 
which modify the lens potentials. Unfortunately, little is actually known about the environments of most 
lenses, so environmental effects in lens models are often unconstrained and sometimes ignored. We show that 
such poor knowledge of environments introduces significant biases and uncertainties into a variety of lensing 
applications. Specifically, we create a mock poor group of 13 galaxies that resembles real groups, generate 
a sample of mock lenses associated with each member galaxy, and then analyze the lenses with standard 
techniques. We find that standard models of 2-image (double) lenses, which neglect environment, grossly 
overestimate both the ellipticity of the lens galaxy (Ae/e ~ 0.5) and the Hubble constant (Ah/h ~ 0.22). 
Standard models of 4-image (quad) lenses, which approximate the environment as a tidal shear, recover the 
ellipticity reasonably well (|Ae/<? < 0.24) but overestimate the Hubble constant (Ah/h ~ 0.15), and have 
significant (^30%) errors in the millilensing analyses used to constrain the amount of substructure in dark 
matter halos. For both doubles and quads, standard models slightly overestimate the velocity dispersion of the 
lens galaxy (Aa/a ~ 0.06), and underestimate the magnifications of the images (A/i//i ~ -0.25). Standard 
analyses that use the statistics of lens populations to place limits on the dark energy overestimate fl\ (by 
0.05-0.14), and underestimate the ratio of quads to doubles (by a factor of 2). The systematic biases related 
to environment help explain some long-standing puzzles (such as the high observed quad/double ratio), but 
aggravate others (such as the low value of Ho inferred from lensing). Most of the biases are caused by neglect 
of the convergence (gravitational focusing) from the mass associated with the environment, but additional 
uncertainty is introduced by neglect of higher-order terms in the lens potential. Fortunately, we show that 
directly observing and modeling lens environments should make it possible to remove the biases and reduce 
the uncertainties. Such sophisticated lensing analyses will require finding the other galaxies that are members 
of the lensing groups, and measuring the group centroids and velocity dispersions, but they should reduce 
systematic effects associated with environments to the few percent level. 

Subject headings: cosmological parameters — dark matter — galaxies: clusters: general — galaxies: halos — 
gravitational lensing 
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1. INTRODUCTION 

The study of strong gravitational lenses offers unique con- 
straints on the masses and properties of galaxy dark matter 
halos (e.g.. [K ochanek l 119911 i Keeton. Kochanek. & Falco 
1998; Rusin, Kochanek, & Keeton 2003b; Treu & Koopmans 
2004), the pr o pertie s of quasars (e.g.7|Nemiroff 1988; 
Richards et al. 200- p) and their host galaxie s (e.g., 
Rixetal.1 I200U Kochane k. Keeton. & McLeod 200 1 ; 
Peng et all r2 004). the Hubble constant (e.g.. iRefsdall 119641 
Kochanek & Schechter 2003) , the nature of dark matter 
(e.g., Metcalf & Madau 2001; Dalai & Kochanek 2002), and 
the properties of dark energy (e.g., Turner 1990; Kochanek 
Il996at IChael 12003: Mitchell et aT] l2004l iLinder 2004). The 
number of strong lenses is approaching 100, with useful 
subsamples that number in the tens (e.g ., Browne et al. 
120031 [Rusin et aTll2003atlO~fek. Rix. & M aoz 2001, and will 
grow dramatically with ongoing and future surveys (e.g., 
Kuhlen. Keeton. & Madau 2004). The relative positions 
and fluxes of the lensed images are routinely measured 
to high precision with th e Hubble Space Telescope and 
radio interferometers (e.g., iLehar et alTl2000t iPatnaik et al.l 
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119991 iTrotter. Winn. & Hewitt! |2000); Einstein rings 
arcs are pro ving to be c ommon in high-resolution near-IR 
images (e.g., Kochanek et al. 2001); and image time delays 
and lens galaxy velocity dispersio ns are succumbing to 
concerted observational effort (e.g., Schechter et al. 199" 
Burudetal. 2002a b; Fassnacht et all 120021 ICollevetal 
20031; iTreu & Koopmans 2002al 120041 With such extensive 
and high-quality data in hand, the results of lensing analyses 
are limited mainly by systematic uncertainties in the lens 
models required to interpret the data. 

We are interested in systematic uncertainties related 
to the fact that lens galaxies are not isolated. Many 
lenses are produced by early-type galaxies that reside 
in overdense regions li ke poor groups o r rich clus- 

ters of galaxi es (e. g., Young et al. 1980; Kundi^yjlJ 

1997a lbl lTonrvlll998t rTbn rv & Kocha n"eklll999t iKneib etaA 
2000; Fassnacht & Lubinl 120021 iMomcheva et al. 2004; 
Williams et al. 2004). The additional mass near the lens con- 
tributes both a "convergence" (additional gravitational fo- 
cusing) and "shear" (gravitational tidal force) 2 to the lens 

2 This shear is formally equivalent to that probed with weak lensing anal- 
yses of clusters and cosmic shear (se e reviews by Bartelmann & Schneider 
1200 II Ivan Waerbeke & Melliei 2003), even though it is detected in a differ- 



potential, plus higher-order terms that may or may not be 
small (see § |2j. One set of systematic uncertainties in 
lens models occurs because the convergence is degenerate 
with the lens galaxy mass, in all lens observables but the 
time delay s; this is the well-known " mass-sheet degener- 
acy" JGorenstein. Shapiro. & Falcolll988t ISahal 12000). Be- 
cause of the degeneracy, the convergence is usually omitted 
from lens models, which leads to a variety of biases. At- 
tempts to correct the biases — for example, by using inde- 
pendent weak lensing analyses to constrain the convergence 
(e.g., Bernstein & Fischer 1999) — have been rare. 

A second set of systematic effects occurs because the 
shear is approximately degenerate with t he elli pticity of the 
lens galaxy (Keeton, Kochanek, & Seliak 1997). Two-image 
(double) lenses suffer badly: even when environmental ef- 
fects are strong, models without shear fit the data just fine, 
and models that include shear reveal it to be degenerate with 
ellipticity. Without independent knowledge of the environ- 
ment, the natural choice for doubles is to omit shear, but 
that must introduce errors into the models. The situation is 
better for four-image (quad) lenses, because the degeneracy 
between ellipticity and shear is only approximate. Both ef- 
fects produce quadrupole terms in the lens potential, but with 
different radial dependences: <f> oc -f r 2 cos 2(9-9 ^) for shear, 
and <j) oc ercos2(9-9 e ) for ellipticity. Quads generally have 
enough constraints to detect the difference, and thus to re- 
veal that shear cannot be neglected (e.g., Keeton et al. 1997). 
Even so, models of quads may not uniquely constrain both the 
ellipticity and shear, leaving a large range of allowed mod- 
els. Furthermore, for both doubles and quads we must ask 
whether approximating the environmental effects as a simple 
shear (neglecting higher-order terms in the lens potential) is 
adequate. 

These issues can be collected into a general question: // 
lens galaxy environments are poorly known, how wrong will 
standard lensing analyses be? We address the question by 
placing galaxies in simple but realistic environments whose 
properties we understand and control, and using them to gen- 
erate catalogs of mock lenses. We then apply standard lensing 
analyses to the mock lenses, and study whether the results ac- 
curately recover the input parameters. We consider a wide 
range of astrophysical and cosmological problems to which 
lensing is applied, including the shapes and masses of galaxy 
dark matter halos, the amount of substructure in dark matter 
halos, the properties of lensed sources, the Hubble constant, 
and the dark energy density. In this paper we begin by exam- 
ining a single (but typical) case to identify problems that arise 
when lens environments are unknown. In a subsequent analy- 
sis we will study in more detail how the errors depend on the 
properties of the environment. 

We stress that we focus on arcsecond-scale lens systems in 
which the lens potential is dominated by one galaxy (o r occa- 
siona l ly two or three galaxie s in close proximity; Rus in et all 
1200 ll iKeeton & WinrJl2003l) . and the contribution from the 
environment can be treated as a perturbation. This is the 
case for nearly all multiply-imaged quasars an d radio sources. 
It differs from the case of lensed arcs (e.g., Gladders et al. 
2003; Zaritskv & Gonzalez 2003) and wide-separation lensed 
quasars (e.g., Inadaetal. 2003), which are produced by the 
large dark matter halos of clusters. 

In this paper we focus on systematic effects associated 
with the angular structure of the potential. There is an ad- 
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ditional set of systematic effects associated with the fact 
that in models of many lenses the radial density profile of 
the lens galaxy is degenerate. We explicitly neglect this 
degeneracy by always assuming that the galaxies can be 
modeled as isothermal ellipsoids. This simple assumption 
seems to be rem arkably good ( e .g.. Izlir itskv & White 1994; 
| Rix et alJ fl997t iGerhardeTail 1200 it McKav eTalT 12002 
Treu & Koopmans 2002a; Koopmans et al. 2003; Rusin et al. 
2003b; Sheldon et al. 2004). Even if this assumption comes 
into question (see Romanowskv et al. 2003), the crucial point 
is that we use a fixed radial density profile for both generating 
and modeling mock lenses. In this way we avoid any system- 
atic errors associated with the radial profile, and highlight sys- 
tematic effects associated with environment. In the follow-up 
analysis we will also examine systematics related to the ra- 
dial profile, and consider whether knowledge of lens galaxy 
environments can actually help break the profile degener- 
acy. Other recent investigations have studied the radial pro- 
file degeneracy in considerable detail (e.g., Treu & Koopmans 
2002a b, 2004; Rusin et al. 2003b). 

This paper is organized as follows. We first create a typical 
lens environment, a poor group of galaxies that mimics the 
group around the observed lens PG 1115+080 (§0. We use 
the galaxies in the group to create a catalog of mock lenses 
(§Ej. We then analyze the mock lenses with standard tech- 
niques to identify environment-induced uncertainties and bi- 
ases in a variety of lensing applications (§ @J. Finally, we 
show how knowledge of the environment can be used to re- 
move the systematic effects (§ [5jl. We offer some general 
comments in §|6j and summarize our conclusions in §0 

2. A MOCK GROUP OF GALAXIES 

Although the full distribution of lens galaxy 
environments is not well kn own, predictions 
( Keeton, Christlein, & Zabludoff 2000b) and observations 
(e.g., Kundic et al. 1997a b; Tonrv 1998; Tonrv & Kochane 
1999; Fassnacht & Lubinl 12002: M omcheva et alJ I20O 
Williams et al. 2004) indicate that the most common lens 
environments are probably poor groups of galaxies. We 
seek to create a mock group that mimics the group at 
redshift z = 0.31 around the 4-image lens PG 11 15+080 
l lKundicetal] fl997a: Tonrv 1998; Momcheva et al. 2004). 
This group appears to be typical of lens environments, and 
to be similar to X-ray se lected grou ps in the nearby universe 
( Zabludoff &Mulchaey 1998, 2000) in terms of its galaxy 
population, kinematics, and X-ray properties (Grant et al. 
2003; Momcheva et al. 2004). In addition, this lens is a 
favori t e for many lensing applicat i ons (e.g., SchechterejjdJ 
[T997t [Kee ton & Kochan ekl Il997[ ISaha & Williams! 119971 
Impevetal. 1998; Zhao & Pronk 2001; Treu & Koopmans 
2002b), so it offers an insightful example. 

The lens in PG 1115+080 consists of four images of a 
quasar at redshif t z s = 1.72 around an ea rly-type galax y at 
redshift z, =0.31. iKundic et"ai] ( ll997al) and lfonrvl ( fl998T) dis- 
covered that there are four other galaxies at the same redshift 
spanning 24" on the sky a nd forming a group wi th velocity 
dispersion a ~ 300 km s" 1 . Momcheva et al. (2004) have ex- 
panded the group membership to 13 galaxies spanning 4.5', 
and calculated a velocity dispersion of a = 354 ±53 km s . 
We use the relative positions of these galaxies, together with 
the relative magnitudes given by Williams et al. ( 2004). 

In deciding how to distribute the mass within the group, 
it is useful to consider two extreme cases that bound the 
possibilities. One case places all of the mass in the group 
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FIG. 1 . — (Left) Logarithmic surface mass density map for the mock group, in the model with all the mass associated with the galaxies and their dark matter 
halos. The map is 6' (570 h kpc) on a side. (Right) Lensing critical curves in the central region of the group. The solid curves show the actual critical curves; 
for comparison, the dotted curves show the critical curve each galaxy would have if it were isolated. The galaxies are labeled with the indices from Tablefll 



member galaxies and their individual dark matter halos; this 
could represent a system where the galaxies have not inter- 
acted tidally with each other. The other case places all of 
the mass in a common dark matter halo, such that the galax- 
ies are essentially massless tracers of the group potential. 
In nearby groups, the kinematics of group member galax- 
ies and the presence of extended X-ray luminous halos sug- 
est that the most of the mass is in a common group halo 
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Zabludoff & Mulchaev 1998; Mulchaey & Zabludoff 1998). 
The situation is less clear for the distant groups (0.3 < z < 1) 
where lens galaxies are likely to be found. X-ray halos 
have been detected in two lensing groups at z ~ 0.3 (includ- 
ing PG 1 1 15+080; Grant et al. 2003), but studying them well 
enough to probe the group mass distributions is exceedingly 
difficult (the X-ray halos are faint extended sources that have 
the bright quasar image superposed). Theoretically, one might 
expect distant groups to be less dynamically evolved than 
nearby groups, and thus to have less mass in a common halo. 
Given the uncertainties, we have analyzed both of these ex- 
treme cases. The systematic errors in lens models are qualita- 
tively quite similar for the two cases, so we present only the 
case where all of the mass is in the galaxies and their individ- 
ual dark matter halos. This choice makes it easier to discuss 
how observations of lens environments can be used to reduce 
or remove the model errors (see §[5}- 

We model the galaxies as isothermal ellipsoids, which 
is consistent with evidence from strong lensin g, stella r 
dynamics, and X-ray studies on small scales JFa bbian' 



1989; Rix et al. 1997; Gerhar d et all2001tlTreu & Koo pman 
2002a; Koopmans et al. 2003; Rusin et al. 2003b), and with 
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[alaxy-galaxy lensing and satellite kinematics on large scales 
Zaritskv & White! [1.994; McKav efal] 120021 iSheldon et all 
2004). An isothermal ellipsoid has projected surface mass 
density (in units of the critical density for lensing) 
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where q < 1 is the axis ratio, 9o is the orientation angle (de- 
fined as a position angle measured East of North), and b is a 



b = 4n 



(2) 



mass parameter related to the velocity dispersion a by 

2 £^ 
D os ' 

where D os and D/ s are angular diameter distances from 
the observer to the source and from the lens to the 
source, respectively. For a spherical galaxy b equals 
the Einstein radius, while for a nonspherical galaxy 
b and R e { n differ by no more than a few percent 
for reasonable axis ratios. The lensing properties 

of iso th ermal ellipsoids are given bv | Kas siola & Kovner 
( 1993), Kormann, Schneider, & Bartelmann ( 1994), and 
Keeton & Kochanek ( 1998). 

The spectra of group member galaxies indicate that many 
are early-type galaxies (Momcheva et al. 2004), so we use the 
Faber-Jackson relation (L oc er 4 ) to estimate the galaxies' ve- 
locity dispersions and b parameters. Specifically, the ratio of 
the b parameters for two galaxies i and j is estimated from 
their relative magnitudes, 



— = lQ-0-2(m,-m } ) _ 

bi 



(3) 



Scatter in the Faber-Jackson relation (e.g. JSheth et alJl2003l) 
could cause scatter between these b ratios and the actual val- 
ues for the PG 1 1 15+080 group, but that is irrelevant because 
we use eq. to define the mock group. We determine the 
scale by setting b = 1 "0 for the galaxy that is an analog of the 
lens galaxy in PG 11 15+080 (galaxy #1 in Table[T|below); this 
value is both consistent with the observed lens (Imp evet alJ 
1998) and typical of all lenses (see the image separation d is- 
tributions in, e.g. JKeeton et all2000bT;lKuhlen et all2 004). 

We assign all of the galaxies an ellipticity e = 1 -q = 0.3, 
which is t ypical for early- type galaxies (e.g ., Bend er et alJ 
ll989USagria et alJll993t [T^rgensen ~alJll995h . and we give 
them random orientations. The resulting parameters for the 
mock group galaxies are given in Tabled A visual sense 
of the group is given in Figure ^ Note that the left panel 
shows the projected mass (not light) distribution, so it is per- 
haps more comparable to maps from A^-body simulations than 
to images of observed groups. 



Table 1 . Properties of the Mock Group 
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NOTE. — Each galaxy has position (x,y), ellipticity and position angle (e, 9q), and mass parameter b\ and feels a convergence K C nv and shear 7 cn v from its 
environment; N 2 and N4 give the number of mock double and quad lenses in our catalog that are associated with each galaxy. 



It is worth pointing out that even though we study a sin- 
gle global system (the group), each galaxy has a unique local 
environment so we actually have 13 sample lens galaxy envi- 
ronments. As a further test of whether the group is reasonable, 
we can quantify the range of environmental contributions to 
the lens potential. For galaxy i, the environmental piece of the 
potential is the part due to all the other galaxies: 



■ = (pi + < 






(4) 



The lensed images tend to lie near a radius (the Einstein ra- 
dius) that is small compared with the typical distance between 
galaxies. Therefore we might expand 4> em j as a Taylor series 
in polar coordinates centered on galaxy i. The lowest order 
significant terms are 3 
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where K e nv and 7 env represent, respectively, the convergence 
(gravitational focusing) and shear (tidal perturbation) dis- 
cussed in the Introduction. Equivalently, the convergence and 
shear can b e expressed in terms of d erivatives of the lens po- 
tential (e.g., Schneider, Ehlers, & Falco 1992): 
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Even if the higher-order terms in eq. (|5) are not negligible, 
the convergence and shear are still a useful way to quantify 
the environmental contribution to the lens potential. 

Table[2lists the convergence and shear computed for each 
galaxy in the mock group, and Figure |5] shows a histogram 
of the shears. For comparison, Holder & Schechter (2003) 
have computed the distribution of shears at the positions of 
early-type galaxies in JV-body and semi-analytic models of 

3 The O(r ) term represents the unobservable zero point of the potential, 
while the 0(r) terms represent an unobservable uniform deflection. 
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FIG. 2. — The histogram shows the shear distribution for the 13 galaxies 
in the mock group. The dashed curve shows the distribution for early-type 
galaxies in galaxy formation simulations 1 Holder & Schechter 2003). 



galaxy formation. This distribution, which is also shown in 
Figure |3 represents the range of shears expected in stan- 
dard galaxy formation models, and is also broadly consistent 
with the shears required to fit observed 4-im age lenses (see 
iKeeton et al.lll997tlHolder & S chechter 2003); so it provides 
at least a rough guide to the shears felt by real lenses. Al- 
though the sample size is small, the shear distribution for our 
mock group is remarkably similar to the distribution found by 
Holder & Schechter. 4 While this test is by no means conclu- 
sive, it does suggest that our mock group is not way off base 
and is in fact quite compatible with the expected distribution 
of environmental contributions to lens potentials. 

3. A MOCK LENS CATALOG 

3.1. Basic approach 

While our mock group represents a single global system, 
the 13 member galaxies all have different local environments 
and therefore sample the range of environments that may be 
common to group galaxies. We have already made use of this 
fact in examining the distribution of shears and using it to ar- 

4 The point in the high-shear tail is galaxy #5, a small galaxy (near the 
center of FigureQ that lies very close to, and is strongly perturbed by, galaxy 
#3. This galaxy has a small lensing cross section, so even if its shear is 
unusually large it contributes little to our analysis and does not affect our 
conclusions. 



gue that our group is reasonable. The next step is to sample 
the range of lens systems associated with these different en- 
vironments. We do this by creating a catalog of mock lenses: 
we place many random sources behind the group, determine 
which are strongly lensed, and tabulate the positions, magni- 
fications, and time delays of all the lensed images. 

We then employ the mock lens catalog to test the effects of 
environment in two ways. First, we imagine that we have ob- 
served the mock lenses, and we fit each one individually with 
standard mass models to constrain the lens galaxy properties 
(mass, ellipticity, substructure, etc.) and Hubble constant Hq. 
Second, we use the statistics of the lens sample to constrain 
the cosmological constant fl\ and to study the relative num- 
bers of quad and double lenses. The main question is whether 
or not standard lensing analyses, in which environment is un- 
known or poorly understood, yield accurate constraints on the 
various parameters. In essence, we are using a Monte Carlo 
approach to sample the errors in lensing analyses that are as- 
sociated with group environments. The remainder of this sec- 
tion is devoted to the details of the Monte Carlo techniques. 

3.2. Application to individual lens systems 

In creating the mock lens catalog, solving the lens equa- 
tion to find the images of a given source is straightforward 
using the algorithm and software by Keeton (2001). The only 
subtleties relate to how the random background sources are 
selected. To sample lens properties realistically, we must ac- 
count for "magnification bias," or the fact that highly magni- 
fied lens systems are more readily discovered than less mag- 
nified systems (e.g., Turner, Ostriker, & Gott 1984). If we 
want to mimic real samples, and in particular if we want to 
avoid undersampling quad lenses, when choosing sources we 
should give more weight to those with large magnifications. 
This amounts to applying a non-uniform probability density 
in the source plane, which has the form 



p(u)du = f(u) 



N(S /fi(u)) 
N(S Q ) 



du, 



(10) 



where /(u) is 1 if a source at position u is multiply-imaged 
and otherwise, /i(u) is the magnification of a source at that 
position, N(S) is the number of background sources brighter 
than flux S, and So is the flux limit of a survey. If the luminos- 
ity function of sources has a power law form, dN/dS oc S~ v , 
then the probability density simplifies to 



p(u)du = f(u) f i(uy- 1 du. 



(11) 



Formally, we just draw from this probability density to gener- 
ate the sample of sources, and then solve the lens equation to 
obtain the mock lens catalog. 

In practice, the probability density eq. d 1 1 1> is quite compli- 
cated (because of the lensing caustics, for example), and there 
is no general algorithm for drawing from a complicated two- 
dimensional probability density. Fortunately, a simplification 
is possible. In the special case i/ = 2we can transform the 
probability density as follows: 



p(u)du = f(u)^2 



;=i 



du 



du = /(x) dx . 



(12) 



where /(x) is 1 if an image at position x corresponds to a 
source that is multiply-imaged, and otherwise. The first 
equality uses the fact that the total magnification is the sum of 
the magnifications for all the images, and the magnification of 



image i is |9x,/9u| (see, e.g., Schneider et al. 1992). The sec- 
ond equality uses the fact that <9x/<9u| is the Jacobian of the 
transformation between the source and image planes. While 
the sum makes the middle expression look complicated, it just 
ensures that the probability density is carried into every part 
of the image plane, with no gaps or overlaps. 5 The bottom line 
is that with v = 2, the magnification weighting in the source 
plane corresponds to uniform weighting in the image plane. 
Therefore, what we can do is sample image positions uni- 
formly within the multiply-imaged region of the image plane, 
and then map them back to the source plane to obtain a set of 
random sources with magnification weighting. 

Apart from ease of use, this transformation has one addi- 
tional advantage and one small drawback. The drawback is 
that we have had to assume a source luminosity function of 
the form dN/dS oc S~ 2 ; but this is not so different from the 
luminosit y fun ction of sources in the largest existing lens sur- 
vey (see § I3.3> . The practical goal of including magnification 
bias is to avoid undersampling quad lenses, and for this pur- 
pose the v = 2 luminosity function is perfectly adequate. The 
fringe benefit of sampling uniformly on the sky is that each 
galaxy in the group is automatically weighted by its lensing 
cross section, so massive galaxies contribute the most to our 
sample of mock lenses just as they would in reality. 

In creating the final catalog, choosing the source density 
represents a compromise between adequate sampling of the 
range of lens properties and the computational time required 
for the lens modeling. We find that a sampling density (in 
the image plane) of 40/D" yields a total of 2606 doubles and 
1043 quads, which seems sufficiently large but still manage- 
able. The number of mock lenses associated with each of 
the 13 galaxies is given in Tabled Once we have created 
the mock lenses, we imagine observing them with the fol- 
lowing measurement uncertainties: 0"003 in the positions, 
which is typical of modern data (from Hubble Space Tele- 
scope images or radio interferometry, see Lehar et al. 2000; 
Patnaik et al. 1999; Trotter et al. 2000); 10% uncertainties in 
the fluxes; and time delay uncertainties of 2 da ys, which is 
somewhat better than much c urrent data (e.g., Schecht er et alJ 
1997; Burudetal. 2002a b) but achievable with dedicated 
campaigns (e.g., Fassnacht et al. 2002; Collev et al. 2003). 

The mock lenses require that we specify the cosmology and 
the redshifts of the source and lens. We adopt a cosmol- 
ogy with Q,m = 0.3, Qa = 0.7, and h = 0.7, and we use the 
lens redshift zi = 0.31 and source redshift z s = 1-72 to mimic 
PG 1115+080. Note, however, that these particular values 
only affect the time delays. The time delay between images i 
and j has the form: 
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The term in square brackets is expressed in angular units, so 
the only physical scale (and the only dependence on redshifts 
and cosmology) appears in the factors out front. Changing 
the redshifts or the cosmology would rescale all time delays 
in the same way, and would not affect our conclusions. 

Once we have the mock lenses, we fit them using the two 
standard lens models that are usually applied when nothing 

5 The way to think about this is that the full source plane maps into the full 
image plane, and each image position corresponds to a single source. 



is known about environment. When a new lens is discov- 
ered, the first thing people usually do is fit it with a singu- 
lar isothermal ellipsoid (SIE). 6 If the fit is poor, the common 
next step is to add an external shear — to model a possible 
environmental perturbation using just the 7 env term in eq. (|5}. 
(The higher-order terms are dropped for simplicity and be- 
cause they are thought to be small; the K env term is dropped 
because it leads to the mass-sheet degeneracy.) Lenses are 
occasionally treated with more complex environment models 
(e.g., Keeton & Kochanek 1997; Bernstein & Fischer 1999; 
Koopmans et al. 2000; Kneib et al. 2000), but those models 
are customized to detailed observational data on the environ- 
ments. Our goal is to identify the kinds of problems that arise 
from standard analyses when little or nothing is known about 
lens environments. 

3.3. Application to lens statistics 

For studying the statistics of lens samples, it is useful to 
make two small modifications to the mock lens catalog. First, 
in this analysis we are not limited by the computational ef- 
fort of modeling the lenses; all we want to do is generate the 
catalog and then determine its statistical properties. We can 
therefore increase the sampling density to ~ 10 /□" to im- 
prove the precision of the Monte Carlo calculations. 

The second modification relates to magnification bias. We 
would like to mimic the largest exi sting lens survey, th e Cos- 
mic Lens All-Sky Survey (CLASS; Browne et al. 2003]). The 
luminosity function of sources in CLASS is wel l described 
as a power law d N/dS oc S~ u with v = 2.\ (Rusin & Tegmark 
l2001l lChae 2003). This is sufficiently simil ar to v = 2 to jus- 
tify the simplifying transformation used in § !3.2l (see eq.ll2>. 
However, since we are modifying the catalog for the statisti- 
cal analysis anyway, we go ahead and use v = 2.1 to allow a 
more direct comparison with CLASS. 7 

The first statistical step is to sum over the sources associated 
with a particular galaxy to obtain that galaxy's lensing cross 
section. In fact, the quantity of interest is the product of the 
lensing cross section and the magnification bias, which we 
term the "biased cross section," and which can be written as 

/v(VMu)) 



du, 



(14) 
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where the integral extends over the multiply-imaged region 
of the source plane, and N(S) is again the number of sources 
brighter than flux S. For a power law source luminosity func- 
tion dN/dS oc S~ y , eq. dl4> simplifies to 
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By integrating over the full multiply-imaged region behind a 
galaxy, we obtain the total biased cross section. With simi- 
lar integrals restricted to the doubly-imaged and quadruply- 
imaged regions, we can compute the biased cross sections for 
double and quad lenses separately. 

The next statistical step is to sum over galaxies to compute 
the lensing optical depth r, which represents the total proba- 
bility that a given source is lensed. In general, we would con- 
sider some appropriate population of galaxies described by 

6 Saha & Williams (1997, 2004) have introduced general non-parametric 
lens models that do not require the isothermal assumption, or indeed any 
other explicit assumptions about the lens galaxy mass distribution. However, 
those models treat the environment as a simple shear, so they are still subject 
to all of the systematic effects that we identify. 

7 Usi ng v = 2.1 prevents us from repeating the coordinate transformation 
used in $ 13.21 so we select points uniformly in the source plan e and apply the 
magnification weighting a posteriori, as indicated by eq. 
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where the first integral is over the volume between the ob- 
server and the source. Replacing a with the biased cross sec- 
tion for doubles or quads, we can compute the optical depths 
T2 and T4 for 2-image and 4-image lenses. The predicted ratio 
of quad to double lenses is then simply Ta/ti. In practice, we 
have only 13 discrete galaxies, so we replace the mass integral 
J(dn/dM)dM with a sum over the group galaxies. Also, we 
know the group redshift, so we drop the volume integral. 

To understand how the optical depth is used to place con- 
straints on Ha, consider a flat universe with a non-evolving 
population of lens galaxies. W hile this simple assumption 
appears to be surprisingly good (Schade et al. 1999; Im et al. 
2002; Ofek et alj l2003t IChae & Mao 2003), allowing some 
evolution would change some quantitative details but not af- 
fect the thrust of our argument (see Mitchell et al. 20QJ). In 
the non-evolving case the optical depth takes the form 



t(z s ) = ax Tn tot x D(z s y 



(17) 
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where w t ot is the integrated comoving number density of 
galaxies, T is a dimensionless number that dep ends on the 
shape of the mass function (see, e.g., Kochanek 1993a; 
Mitchell et al. 2004), and D(z) is the comoving distance to 
redshift z. In other words, the optical depth is proportional to 
the volume of the universe out to the redshift of the sources. 
Conceptually, one takes a measurement of the optical depth, 
adds data or models for the mass function and a model for 
the cross section, and uses eq. dl7> to place constraints on the 
volume of the universe and hence £!a. 

4. IDENTIFYING THE PROBLEMS 
4. 1 . Basic results: \ 2 , e > an d h 

We first consider some of the basic quantities derived from 
models of individual lenses: the goodness of fit \ 2 , the mass 
ellipticity e, and the Hubble parameter h. Figure[3]shows his- 
tograms of these quantities for SIE and SIE+s hear models of 
the mock double lenses in the catalog from § 13.21 Doubles 
provide just 7 observables: position and flux for each image, 
and the time delay. SIE models have 7 parameters: b, e, and 
#o for the galaxy, position and flux for the source, and h. SIE 
models of doubles therefore have A'dof = 0, and in general it 
is possible to find models that fit the data perfectly (% 2 = 0). 
Unfortunately, there are large errors in the recovered mass el- 
lipticity and Hubble parameter caused by the neglect of envi- 
ronment. The problem is worse with SIE+shear models; even 
though environment is included (in the shear approximation), 
it is completely unconstrained and so there are enormous un- 
certainties in the models. Doubles suffer a fundamental prob- 
lem: with so few constraints, lens models cannot even detect 
an environmental component in the lens potential, much less 
constrain it well enough to yield accurate results. 

Figure|4]shows that the situation is better with quad lenses 
because the additional images provide more constraints. 
Quads offer 15 observables: position and flux for each im- 
age, and three independent time delays. So both SIE models 
(7 parameters) and SIE+shear models (9 parameters) are well 
constrained. In fact, SIE models generally give terrible fits: 
among our mock quads the median \ 2 is 464, and in 95% of 
the lenses the model can be ruled out at more than 99% con- 
fidence (x 2 > 20.1 for A^of = 8). In other words, quads have 
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FIG. 3. — Results for models of mock double lenses. SIE models are 
shown in the left column, and SIE+shear models in the right. The panels show 
histograms of \ 2 values (top), the inferred lens galaxy ellipticity (middle), 
and the Hubble constant (bottom). (Note that all values outside the x-axis 
range are placed in the left- and right-most bins.) The input values e = 0.3 
and h = 0.7 are indicated by dashed lines. 



FIG. 4. — Similar to Figure|5] but for mock quad lenses. The open his- 
tograms show all quads, while the shaded histograms show only those where 
the lens model provides an acceptable fit to the data (models that are not ruled 
out at the 99% confidence level; see text). 



enough constraints to make it obvious that environment can- 
not be neglected. In this case the significant errors in e and h 
are irrelevant because the % 2 already reveals that the models 
are wrong. Like real lenses (Keeton et al. 1997), our mock 
lenses indicate that quads cannot be well fit by models that 
neglect environment. 

The SIE+shear models of quads are more interesting. Now 
the median % 2 among the mock quads is 19.3, and in only 
52% of the systems can the model be excluded at 99% confi- 
dence (x 2 > 16.8forAW = 6). In the other 48% the SIE+shear 
model provides an acceptable fit to the data. Thus, approx- 
imating the environment as a shear is sufficient to provide 
a good fit to many quads. This approximation is not ade- 
quate, however, for recovering accurate values for the ellip- 
ticity and Hubble parameter. The errors in these quantities, 
while smaller than for doubles, are still worrisome. 

Particularly troubling is the fact that the errors are not 
random but biased. For doubles, SIE models usually over- 
estimate both e and h: the median recovered ellipticity is 
e = 0.47; and (99,66,30,7)% of the mock doubles yield h > 
(0.7,0.8,0.9,1.0). For quads, SIE+shear models do reason- 
ably well with ellipticity: 88% of the quads with acceptable 
fits (and 74% of all quads) have ellipticity errors |Ae| < 0.1. 
But h is still a problem, as (98,26, 14,7)% of the quads with 
acceptable fits yield h > (0.7,0.8,0.9, 1.0). Poor knowledge 
of environment can clearly cause not just uncertainties but 
also significant biases in lensing measurements of ellipticity 
andHo. 

In the remainder of the paper we use the SIE lens models for 
double lenses (because they give perfect fits, and SIE+shear 
models are underconstrained), and the SIE+shear models for 
quads. 

4.2. Mass and velocity dispersion of the lens galaxy 



It is often said that lensing robustly measures the mass of 
the lens galaxy within the Einstein radius. To be precise, this 
statement applies to quad and Einstein ring lenses, while for 
doubles lensing measures a mass-radius relation of the form 
M(Ri)/Ri + M(R 2 )/R 2 = irE cdt (Ri+R 2 ) where /? , and R 2 are 
the ima ge radii and M(R) is the aperture mass (see Rusin et al. 
2003b for a detailed discussion). (Quads and rings can actu- 
ally be thought of in the same way, but with R\ = R 2 = R e i D 
since all the images appear near the Einstein radius.) For sim- 
plicity we limit our discussion to the mass within the Einstein 
radius as appropriate for quads and rings, keeping in mind that 
it would be modified slightly for doubles. 

An important caveat is that lensing measures only the total 
projected mass within the Einstein radius. The total certainly 
includes contributions from the lens galaxy (its stellar com- 
ponent and dark matter halo). But there is also a contribution 
of 7r^ 2 in K env £ cr jt from the convergence. Physically, this mass 
might represent a smooth background in which the galaxy is 
embedded, such as a common group dark matter halo; or it 
might represent mass in the immediate foreground or back- 
ground, such as the halos of other group member galaxies that 
overlap the line of sight. Models that neglect convergence as- 
sume that all of the measured mass must come from the lens 
galaxy, and hence overestimate the galaxy mass. 

To examine this effect, it is more instructive to quote the ve- 
locity dispersion of the lens galaxy; while the mass depends 
strongly on the aperture and is not directly observable, the ve- 
locity dispersion depends only weakly on aperture (not at all 
for isothermal models) and is observable (actually the veloc- 
ity dispersion of the stars; see below). For isothermal models, 
the mass within the Einstein radius is M oc bR e i n oc a 2 R e i n (see 
eq. |2}, so the conversion is simple. Figure [5] shows that the 
models do slightly overestimate the velocity dispersion, with 
a typical error of ^6% (and a broad range, especially for dou- 
bles). While small, the errors ar e larger than the ty pical ~1% 
measurement uncertainties (see Rusin et al. 2003b), and they 
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FIG. 5. — Errors in the velocity dispersion of the lens galaxy; <7tme is the 
true (input) value, while <T mo d is the value inferred from lens models. The left 
panel shows doubles; the right panel shows quads, where the open histogram 
show all quads while the shaded histograms show only those where the model 
fits the data. The dashed lines shows where the ratio is unity. 



represent a systematic shift rather than random errors. 

The errors may bear on a historic controversy in lens- 
ing analyses: the relation between the velocity dispersion of 
the dark matter (t he a p arameter in isothermal models) and 
that of the stars. iGottl 11977) argued that the ratio should 
be <7DM/cstars = (3/2) 1 / 2 if galaxy luminosity densities are 
r" 3 power laws, but studies of more realistic models for the 
luminosity indicated that the central stellar velocity d isper- 
sion err) actu ally obeys ctdm/co ~ 1 JFranxll 1993t iKochanekl 
1993b, 1994). The difference may seem small in the veloc- 
ity dispersion, but it looms large in lens statistics because 
the lensing optical depth scales asr« 0"q M . Most recently, 
Treu & Koopmans (2004) have found (cdm/co) = 1.15 ±0.05 
by combining direct measurements of <jq for five lens galax- 
ies with model determinations of ctdm- We cannot measure 
this ratio in our simulations since we do not include a stellar 
component. However, we can say that if standard models are 
overestimating <7dm then this effect could explain part but not 
all of the difference of the measured value from unity. 

4.3. Image magnifications 

In some applications lensing is used as a natural tele- 
scope that provides high resolution for studyin g the struc- 
ture of quasars and their host galaxies (e.g., Rix et al. 2001; 
Kochaneket al. 200ir.lPen g et al.l20 04: Richards et al. 2004). 
In this case accurate lens models are essential for "de-lensing" 
the system, or removing the effects of lensing magnification 
and distortion to infer the intrinsic properties of the source. In 
other applications the need to know the magnificatio n an d dis- 
tortion is less obvious but no less important (see § I4.4> . We 
consider how reliable lens models are for such applications, 
focusing on the magnifications of unresolved images since we 
use only point-like sources. 

When studying magnifications it is important to separate 
images of different types: images that lie at minima of the 
time delay surface and have positive parity, as opposed to im- 
ages that lie at saddlepoints of the time delay surface and are 
parity-reversed (see, e.g., Schneid er et alj fl992). A double 
has one minimum and one saddlepoint, while a quad has two 
of each. The image type has become particularl y im portant 
in studies of microlensing and millilensing (see § I4.4> . If the 
lensed images are point-like then the parities cannot be ob- 
served directly, but image parities can almost always be de- 
termined unambiguously with even simple lens models (e.g., 
Saha& Williams 2003). 

Figures [5] and \J\ show histograms of the error in the lens- 
ing magnification, for doubles and quads respectively. We see 
that standard lens models tend to underestimate the magni- 
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FIG. 6. — Histograms of the ratio of the lensing magnification inferred 
from lens models to the true magnification, for mock double lenses. The 
panels refer to the two image classes (saddlepoints and minima). The dashed 
lines show where the ratio is unity. 




FIG. 7. — Similar to Figure|6| but for mock quad lenses. The panels refer to 
the four image classes. The open histograms show all quads, while the shaded 
histograms show only those where the lens model provides an acceptable fit 
to the data. 



fications, often by a factor of ~1.5—2 and sometimes much 
more. The important qualitative point is that the errors are 
biased: the models almost always underestimate the magnifi- 
cations. 

4.4. Constraints on small-scale structure 

Lensing offers unique constraints on small-scale struc- 
ture in lens galaxies via the ph enomena of microlens- 
ing by stars (e.g., |C hang & Refsd alj |1979t iPaczvnskil 
1986; Schmidt & Wambsganssl Il998t iWvithe et alJ 
2000) and millilensing by dark matter 
Mao & Schneider 1998; Metcalf & Madau 2001; 
2002; Dalai & Kochanek 2002). The question is whether 
environment-related errors in lens models can affect the 
results. The millilensing constraints, in particular, are derived 
from lenses with "flux ratio anomalies," or flux ratios that 
are inconsistent with smooth lens potentials. Conceptually, 
the probability of having a flux ratio anomaly is given by the 
abundance of mass clumps times the cross section for any 
given clump to cause millilensing (to significantly perturb 
the flux of an image). Converting the abundance of flux 
ratio anomalies to a clump abundance therefore requires 
the ability to compute a clump's millilensing cross section. 
But that cross section depends in a complicated way on the 
total lens potential. The dependence is not apparent in the 
detailed simu l ations used by [Dalai & Kochanek (2002) and 
Metcalf et al. (2003) to derive clump abundances, but it is 
made explicit in the analytic expressions derived by Keeton 
(2003) for the millilensing cross sections of clumps modeled 
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FIG. 8. — Errors in the estimated cross section for millilensing by a mass 
clump in front of an image in a mock quad lens. The histograms show the ra- 
tio of the cross section computed with the fitted lens model, to that computed 
with the true lens potential. The dashed lines show where the ratio is unity. 
The open histograms show all quads, while the shaded histograms show only 
those where the lens model provides an acceptable fit to the data. 



as isothermal spheres. 

We focus on millilensing because there are as yet no sim- 
ple expressions for microlensing cross sections. We compute 
the cross section for an isothermal clump to produce a 30% 
change in the flux of an image. Figure [8] shows the ratio of 
the cross section computed with an SIE+shear lens model to 
the cross section computed with the true lens potential, for the 
mock quad lenses. 8 Studying the four types of images sepa- 
rately is important because saddlepoint images — especially 
the bright saddles — are believed to be particularly suscep- 
tible to perturbations; this effect provides a unique signature 
of lensing by small-scale stru cture that allows it to be distin- 
guished from other effects (Schechter& Wambsaanss 2002; 
Kochanek& Dalai 2003). 

There are indeed errors in the millilensing cross sections 
due to errors in the lens models. Interestingly, the error distri- 
butions appear to be symmetric in ln(A mo d/A true ). They tend 
to have relatively tight cores but broad wings, so the RMS er- 
ror in ln(A mo d/A tlue ) is 0.29 for the bright saddles, 0.31 for the 
faint saddles, 0. 14 for the bright minima, and 0.09 for the faint 
minima. (We are considering only quads where the SIE+shear 
model provides an acceptable fit.) The fact that bright saddle- 
point images are so sensitive to the model errors is troubling 
because they are so important for millilensing analyses. 

One would naively expect that a ^30% uncertainty in the 
millilensing cross section might cause a comparable uncer- 
tainty in the inferred clump abundance, which would make 
this effect smaller than other statistical and systematic uncer- 
tainties in current results (see Dalai & Kochanek 2002). How- 
ever, to make that statement quantitative and precise we would 
need to replicate complete millilensing analyses (following, 
e.g., Dalai & Kochanek 2002; Metcalf et al. 200J). As for 
microlensing, without a simple way to compute perturbation 
cross sections we can only guess that there would be qualita- 
tively similar results. Clearly there is much more to do here; 
we just want to point out that analyses of small-scale structure 
depend on the overall lens model, and errors in that model re- 

8 Small-scale structure analyses are usually restricted to quads, because 
doubles have too few constraints to allow robust identification of flux ratio 
anomalies. 




FIG. 9. — Histogram of errors in the lensing cross section; cr mo & is the 
biased cross section with the galaxy assumed to be isolated, while <T truc is the 
value with the galaxy in its proper environment. 



lated to environment cannot be ignored. 

4.5. Constraints on FIa 

We now switch gears and turn to statistical analyses of lens 
populations and the constraints they yie ld on the cosmologi- 
cal constant flj\. As discussed in § 13.31 (see eq.U7K lensing 
limits on fl\ depend on the optical depth t, so the key ques- 
tion is whether environment causes errors in r. First, though, 
it is useful to ask whether there are errors in the lensing cross 
section (from which the optical depth is determined). To an- 
swer this question, we compare the biased cross section a^ue 
that a galaxy has when placed in its proper environment to 
the biased cross section <r mo d the same galaxy would have if 
it were assumed to be isolated. Any differences indicate that 
neglecting environment does cause errors in the cross section. 
Figure [9] shows the errors in the form of a histogram of the 
ratio (7 mo( j/tTtrue. Models that neglect environment underesti- 
mate the cross section by tens of percent or more. The un- 
derestimate is caused entirely by the magnification bias, 9 and 
arises beca use s tandard models underestimate the magnifica- 
tions (see § I4.3> . 

The next question is how errors in the cross section propa- 
gate into errors in 17a- Figure [H)| gives a sense of the effect, 
by showing the error in tt\ that would result from each of the 
cross section errors in Figure [9] Underestimates in the cross 
sections lead directly to overestimates in £Ia, which can be 
quite large. This figure is a bit unfair, because constraints on 
fl\ properly depend on the optical depth, and not all galaxies 
contribute equally. For example, the point at fl\ ss 1 (corre- 
sponding to the point with CT mo d/<T tl -u e = 0.36 in Figure |9j is 
produced by galaxy #5, which is strongly perturbed by galaxy 
#3. However, galaxy #5 has a small cross section and con- 
tributes little to any lens sample, so its properties are not very 
important. A better approach is to compute the optical depth 
by summing the cross sections, and then compare the "model" 
(galaxies assumed to be isolated) and true optical depths. This 
lead to an error r mo d/T tnie = 0.70, and hence an inferred value 
Oa = 0.84 (indicated by an arrow in Figure fTol l. 

9 Curiously, the radial caustic of an isothermal lens is not affected by con- 
vergence and shear from the environment; the radial caustic is the set of points 
that map to the origin, and the deflection from convergence and shear van- 
ishes at the origin (see eq. |5}. Except in the rare case of naked cusps, the 
unbiased cross section is simply the area enclosed by the radial caustic, so it 
is insensitive to convergence and shear. 
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FIG. 10. — Errors in the inferred value of f^A, assuming a typical source 
redshift z s ~ 1.27 as appropriate for CLASS i Marlow et al. 2000; Chae 
2003). The histogram shows the values for the 13 galaxies taken separately. 
The dashed line shows the input value Qa = 0.7. The arrow at Qa = 0.84 
shows the value if we average over the 13 galaxies (weighting each by its 
cross section). The arrow at Qa = 0.75 shows the value if we assume that 
25% of all lens galaxies lie in groups like our mock group, while the other 
75% are isolated. 
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galaxy index 

FIG. 1 1 . — Quad/double ratios for the 13 galaxies, computed as the ratio 
of the biased cross section for 4-image lenses to that for 2-image lenses. The 
stars show the results with each galaxy assumed to be isolated; the galaxies 
should have the same quad/double ratio since they have the same ellipticity, 
so the fluctuations indicate the level of numerical noise (~3%). The squares 
show the results with each galaxy placed in its proper environment. (The 
lines are drawn just to guide the eye.) The arrow at Q/D = 0.50 shows the net 
ratio when all the galaxy cross sections are summed. 



This result holds if all lens galaxies lie in environments like 
our mock group, which is unlikely. The final effect of en- 
vironment on lens statistics cannot be determined because the 
distribution of lens galaxy environments is unknown. Still, we 
can make a crude estimate. Keeton et al. (2000b) use galaxy 
demographics to predict that at least ^25% of lens galaxies 
lie in groups or clusters. If we conservatively assume that 
the fraction of lenses with groups is indeed 25% while the 
other 75% are isolated, we would find a cross section error 
7"mod/Ttrue = 0.90 and hence an inferred value Qa = 0.75 (also 
indicated in Figure 1101 . Even this simple and conservative 
estimate points to the importance of measuring lens galaxy 
environments and understanding their effects on lens statis- 
tics. 

4.6. The quad/double ratio 

A long-standing puzzle in lens statistics is why so many 
lense s are quads rather than do ubles (Kinaetal. 1996; 
Kochanek 1996b; Keeton et al. 1997; Rusin & Tegmark 2001; 
Cohn & Kochanek 2003). In the CLASS statistical sample of 
13 lenses, 7 are doubles and 5 are quads (Browne et al. 2003). 
(The remaining lens is a complicated system that has six im- 
ages because there are three lens galaxies; see Rusin et al. 
2001.) The quad/double ratio depends on the relative cross 
sections for 4-image and 2-image lenses (along with mag- 
nification bias), and it is thought to depend mainly on the 
ellipticity of lens galaxies. An observed quad/double ratio 
near unity is said to require a typical ellipticity e ~ 0.6 that 
is much larger than obse rved or pr edicted fo r normal galax- 
ies (Kochanek 1996b; Keeton et al. 1997; Rusin & Tegmark 



2001). In fact, some have spoken of the high o bserved ra- 
tio cre ating an "ellipticity crisis" in lensing (see Kochanek 

1996B- 

Conventional wisdom holds that environment has little ef- 
fect on the quad/double ratio, because (when averaged over 
orientation) shear does not significantl y change the 2-imag e 
and 4-ima ge cross sections ( e .g.. | Rusin & Tegmarkl 12001). 
However, Cohn & Kochanek (2003) recently showed that 
low-mass satellites around lens galaxies, which are common 
and not very sensitive to the larger environment, can roughly 
double the predicted quad/double ratio. We can now deter- 



mine the effects of a larger group environment, using a model 
that is more realistic than simple shear. 

Under the standard assumption that galaxies are isolated, 
we would compute a quad/double ratio of 0.25 for each of 
our galaxies (the same because they all have the same ellip- 
ticity). By contrast, placing the galaxies in their proper envi- 
ronments yields the quad/double ratios shown in Figure [^ 
We find that neglecting environment causes significant er- 
rors in the quad/double ratio. More relevant than the indi- 
vidual ratios is the net ratio after summing the 2-image and 
4-image cross sections for all the galaxies; this represents the 
ratio that would be expected for a real lens survey. The net 
quad/double ratio is again 0.25 when environment is ignored, 
but 0.50 when environment is included. In other words, ne- 
glecting environment causes a significant underestimate of the 
quad/double ratio. 

Note that our predicted quad/double ratio of 0.50 is higher 
than the ratio 1043/2606 = 0.40 in our mock lenses. Part 
of the difference may be due to statistical fluctuations in the 
mock lens catalog, because some of the galaxies produce a 
small number of lenses (see Table QJ. However, most of the 
difference arises from magnification bias. The mock lenses 
were generated using magnification bias corresponding to a 
source luminosity function dN/dS oc S~ y with v = 2 (see 
§|3jl. By contrast, the lens statistics calculations use a slightly 
steeper source luminosity function with v = -2. 1 (as appropri- 
ate for comparison to CLASS). The steeper luminosity func- 
tion produces a larger magnification bias that increases the 
weight of quads relative to doubles. 

Even with environment included our predicted quad/double 
ratio is still somewhat smaller than observed. Part of the ex- 
planation may be that we have assumed a particular ellipticity, 
e = 0.3. Increasing the ellipticity to e = 0.4 but holding every- 
thing else fixed would yield a quad/double ratio of 0.67 with 
environment included, which is close to the observed ratio. 
(With e = 0.4 and environment omitted, the quad/double ra- 
tio would be 0.40.) A proper analysis would have to include 
an appropriate distribution of ellipticities. Another part of the 
explanation may be that we have considered just one mock 
group, and the full distribution of lens environments will need 
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to be taken into account. A third issue is that we have not con- 
sidered low-mass satellite gal axies, which would furth er in- 
crease the quad/double ratio (Cohn & Kochanek 2003). Still 
a fourth possibility is that the high observed quad/double ra- 
tio is a statistical fluke. Thorough study of these possibilities 
is beyond the scope of this paper. We mainly want to point 
out that environment represents a systematic effect that can 
significantly increase the predicted quad/double ratio without 
violating other constraints (c.f. Rusin & Tegmark 2001). 

5. FIXING THE PROBLEMS 

Having identified some of the environment-related prob- 
lems in lensing analyses, we briefly consider what causes 
them and how they can be fixed. (A fuller treatment of these 
issues will be given in the follow-up paper.) For doubles the 
situation is clear: with so few constraints, lens models are un- 
able to recognize and constrain an environmental contribution 
to the lens potential; to fix the errors, it will be necessary to 
study the environment separately and impose it on the lens 
models. 

For quads, lens models are often able to constrain shear 
from the environment, but the shear approximation clearly 
fails to capture all of the environmental effects. Returning to 
eq. Q, the question is whether the errors are caused mainly 
by neglect of the convergence term (K e „v, representing the ad- 
ditional mass at the position of the lens galaxy), or by neglect 
of the higher-order terms. To address this question, we con- 
sider what would happen if we somehow knew K en v for each 
galaxy and could include it in the models . According to th e 
mass-sheet degeneracy (Gorensteinet al. 1988; Saha 2000), 
the K env term simply rescales some of the model quantities 
without affecting the goodness of fit; some of the key rescal- 
ings are: 



boc(l-K env ), 

croc(l-K env ) I//2 , 
/!0C(1-K e „v), 
^0C(1-K env )" 2 , 
aaa-Kenv)- 2 ^ 1 ', 



(18) 
(19) 
(20) 

(21) 

(22) 

where a is the inferred velocity dispersion while a is the bi- 
ased cross section, and v is the power law index of the source 
number counts (see §0. Figures IT2l and H31 show that this 
rescaling makes the model results much more accurate. In 
particular, it removes the biases in quantities such as the el- 
lipticity, Hubble parameter, lens galaxy velocity dispersion, 
image magnifications, and fl\. There are still some random 
errors: 0.07 in e and h; and 0.20 and 0.16 in ln/i (for the 
bright saddle and bright minimum images, respectively). 10 
Interestingly, adding the convergence does not fix the under- 
estimated quad/double ratio, because it rescales the 2-image 
and 4-image cross sections in the same way leaving the ratio 
unchanged. Still, it appears that neglecting the convergence 
from the environment is the primary cause of the biases in 
lensing results for quads. Whether the residual scatter is due 
to neglect of the higher-order terms in eq. Q or just to obser- 
vational noise is not clear from this analysis. 

The problem is that the convergence is not directly observ- 
able, because it usually represents dark matter. Nor is there 
a reliable way to estimate the convergence from properties of 
lens models. This statement may seem surprising, because 

10 These are the numbers for quads in whi ch th e lens models provide an 
acceptable fit (the shaded histograms in Figure [T2*l . 
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FIG. 12. — Sample results for mock quad lenses, for SIE+shear models 
that also include the convergence K cnv from the environment. We show his- 
tograms of the goodness of fit x 2 > the ellipticity, and the Hubble parameter 
(left), as well as errors in the lens galaxy velocity dispersion and the magni- 
fications of the bright saddle and bright minimum images (right); this figure 
is to be compared with Figures[4]|5] and!71(but note the different axis scales). 
As before, the open histograms show all quads while the shaded histograms 
show those where the model provides an acceptable fit. 
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FIG. 13. — Similar to Figure flOl but for SIE+shear models that also include 
the convergence re cn v from the environment. The dashed line shows the input 
value Qa = 0.7. The arrow at Qa = 0.712 shows the value if we average 
over the 13 galaxies (weighting each by its cross section). The arrow at f2yv = 
0.703 shows the value if we assume that 25% of all lens galaxies lie in groups 
like our mock group, while the other 75% are isolated. 



for simple environments dominated by a single spherical halo 
there is a general relation between shear — which can be in- 
ferred from lens models, at least for q uads — and convergence 
jMiralda-Escudel 1 99 U iKaiserl 1 995t) : 



la 



v)r- 



(23) 



where R is the offset between the lens galaxy and the environ- 
ment halo, and (K e nv)j? is the mean convergence (or surface 
mass density in units of the critical density for lensing) in an 
aperture of radius R centered on the environment halo. Be- 
cause of the average, n em and 7 env depend on radius and on the 
halo properties (mass, concentration, etc.) in different ways, 
and there is no universal relation between them. A second 
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FIG. 14. — Similar to Figure [T2l but for models that incorporate all of the 
group member galaxies using reasonable observational constraints. Here all 
of the quads are well fit by the models. 
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FIG. 15. — Similar to Figure [l4l but for mock double lenses. 



and more troubling problem is that the convergence is a scalar 
while the shear is a rank-2 traceless tensor (or a headless vec- 
tor with an amplitude and a direction that is invariant under 
180° rotations; see eq.[3Jl. This means that in an environment 
comprising multiple mass components, the convergence and 
shear from the different components sum in different ways, 
and attempts to relate the net convergence to the net shear 
break down. A final problem is that this could never work for 
doubles anyway, since there are no good model constraints on 
the shear. 

The only solution is to build a model that allows us to trans- 
late observable quantities into reliable inferences about the 
convergence. If we can identify and model all of the mass 
components in the lens environment, then the shear, con- 
vergence, and higher-order terms will all be included self- 
consistently. Creating such a model is not as crazy as it 
may sound, because many of the key parameters are observ- 
able. For example, the relative positions of the group member 
galaxies can be measured, and the relative magnitudes can 
be used to derive the b ratios (recall the Faber-Jackson rela- 
tion, eq.[3}. Although the mass ellipticities cannot be directly 
observed, we can hypothesize that they are not critical; we 
can set them to zero when modeling the mock lenses and see 
whether we still obtain accurate model results. (To be pre- 
cise, we let the main lens galaxy be elliptical but make all the 
other galaxies circular.) In other words, all of the new pa- 
rameters in the most basic realistic model of environment are 
either constrained by observations or fixed, so adding com- 
plexity has not caused an explosion of free parameters. We 
now consider whether such models are sufficient to fix the 
problems in lensing analyses. To allow for realistic noise in 
the observational constraints on the new parameters, we in- 
clude generous uncertainties: measurement errors of 1" in the 
relative positions; and intrinsic scatter of 0.08 dex in the <r(L) 
relation (Sheth et al. 2003), or 0.16 dex in the model b ratios. 

Figure^Jshows the results when we apply such models to 
quad lenses. There are no biases in the results, and the scatter 



is encouragingly small: just 0.02 in both e and h; and 0.05- 
0.06 in ln/i. In other words, models that include the environ- 
ment break the mass-sheet degeneracy and fix the problems 
with standard lensing analyses. Furthermore, these models do 
notably better than the models in Figures [T2l and [T3l that just 
included convergence, indicating that the higher-order terms 
in the lens potential from the environment (eq. [5} are in fact 
significant. Note that neither the assumption of circular group 
galaxies nor the inclusion of noise causes significant problems 
in the models. 

Figure [TBI shows that similar results hold for doubles, with 
slightly larger scatter: 0.07 in e; 0.03 in h; and 0.08 in ln fi. 
While we saw in §|4]that doubles are largely useless for astro- 
physical measurements when environments are unknown, we 
see now that observing and modeling the environment makes 
doubles almost as valuable as quads. Finally, Figures [TBI and 
[T71 show that fl\ and the quad/double ratio can also be re- 
covered well from statistical models that account for environ- 
ment. 

6. DISCUSSION 

In this paper we have effectively studied 13 different 
galaxy/environment configurations, but they were all drawn 
from a single mock group of galaxies. We must therefore 
acknowledge the following caveats and questions about our 
results. 

Is our mock group realistic? We constructed the sys- 
tem to mimic the group around the observed 4-image lens 
PG 1115+080, and to be similar to nearby X-ray luminous 
groups, so it has reasonable properties. Furthermore, the dis- 
tribution of lensing shears produced by the system is simi- 
lar to what is measured a nd predicted for (quad) lenses (c.f. 
Holder & Schechter 2003), so the group's contribution to the 
lens potential seems realistic. However, the proper way to an- 
swer this question is to study more real groups around real 
lenses in detail. 

Is this group typical? Our group represents just one sample 
system. Based on current knowledge of lens environments 
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inferred value of H A 

FIG. 16. — Similar to Figure fT3l but for models that incorporate the group 
member galaxies using reasonable observational constraints. The arrow at 
SI A =0.720 shows the value if we average over the 1 3 galaxies (weighting 
each by its cross section). The arrow at £7a = 0.705 shows the value if we 
assume that 25% of all lens galaxies lie in groups like our mock group, while 
the other 75% are isolated. 




galaxy index 

FIG. 17. — Similar to Figure llll but for models that incorporate the group 
member galaxies. The squares again show the quad/double ratios when the 
galaxies are placed in their true environments, while the stars show the values 
using the model environments. The net ratio with the model environments is 
0.45 compared with the true value of 0.50. 



(limited though it is) and on the shear distribution, we believe 
that it captures features that are representative. But we must 
identify the full distribution of lens environments before we 
can say for sure. 

How do the results depend on the group properties? It 
seems likely that the systematic effects will scale with the 
mass or richness of the environment. However, that glib gen- 
eralization obscures a lot of interesting details. A group's con- 
tribution to the lens potential may depend on many properties 
beyond its total mass: the ellipticities of the group member 
galaxies (addressed briefly in this paper); the properties of the 
galaxy population, such as the elliptical/spiral and dwarf/giant 
ratios; the fraction of mass in the common halo versus that 
bound to the galaxies (e.g., the degree to which the galaxy 
halos have been truncated by tidal stripping); the radial pro- 
file and angular structure of the common halo; and the offset 
between the group centroid and the lens galaxy. In the limit 
that the environment is dominated by a common dark halo, we 
must also understand how well observables such as the group 
centroid and velocity dispersion will be able to constrain the 
important properties of the environment. Studying these is- 
sues in detail is the subject of the follow-up paper. 

Much more detailed understanding of lens environments is 



clearly required. We need to characterize the distribution of 
lens environments in order to guide theoretical calculations. 
For lenses lying in groups (or clusters), we must find the mem- 
ber galaxies and measure the centroid and velocity dispersion 
of the system, and then build sophisticated lens models that fi- 
nally treat environments properly. Even for lenses that do not 
lie in groups or clusters, careful observations will be needed 
in order to establish that the environments are unimportant. 

7. CONCLUSIONS 

Poor groups of galaxies around strong gravitational lens 
systems can affect lensing constraints on the masses and 
shapes of galaxy dark matter halos, the amount of substruc- 
ture in dark matter halos, the quad/double ratio, the properties 
of lensed sources, the Hubble constant, and the dark energy 
density. Not knowing that a lens lies in a group can cause bi- 
ases and uncertainties in lensing analyses. Models of double 
lenses that neglect the environment will generally fit the data 
well (\ 2 ~ 0) but yield parameter values that are grossly in- 
correct. Models of quad lenses will reveal that environment 
cannot be ignored but will not fully constrain the environ- 
mental component of the lens potential. The standard shear 
approximation for quads is generally adequate for fitting the 
data, but not for recovering correct parameter values. The 
essential point is that the environment-related errors are not 
random errors but systematic biases, such as overestimates 
of lens galaxy velocity dispersions, the Hubble constant, and 
17a, and underestimates of the image magnifications. 

These systematic effects help resolve one long-standing 
puzzle in lensing, but seem to aggravate two others. The 
solved problem is the high number of quad lenses relative to 
doubles (also see Cohn & Kochanek 2003). We have shown 
that neglecting environment can cause models to significantly 
underestimate the quad/double ratio; accounting for environ- 
ment can increase the ratio to near the observed value. 

The first unsolved problem involves lensing and the Hub- 
ble constant. Analyses of lens time delays that invoke what 
we think we know about galaxy dark matter halos yield Hq ~ 
50 km s" 1 Mpc" 1 , which is low compared with the conven- 
tional value Hq ~ 70 (see Kochanek 2002, 2003, and refer- 
ences therein). We find that that poor knowledge of environ- 
ments causes an overestimate of Ho, so accounting for envi- 
ronment would worsen the problem. The extent of this prob- 
lem is actually uncle ar, because some other lensing analyses 
yield H ~ 70 (e.g., Koopmans et al. 2003; Saha & Williams 
2004). In any case, there is a discrepancy remaining to be 
resolved, and the role of environments cannot be ignored. 

The second puzzle involves lensing constraints on fl\. Tra- 
ditional lensing analyses have tended to produce low values 
of Ha, such as the oft-quoted bound Ha < 0.66 at 95% con- 
fidence (Kochanek 1996a), that are in marginal conflict with 
the concordance cosmology. We find that poor knowledge 
of environments leads to an overestimate of Ha, which would 
seem to worsen the problem. However, updated analyses have 
revised the lensing results upwards, with best-fit values in the 
range H A ~ 0.72-0.78 (Mitchell et al. 2004). In this case, cor- 
recting for environment might even resolve any discrepancies 
that may remain between lensing and other methods. 

Apart from a general desire to make lensing analyse s accu- 
rate, there is new impetus to control the systematics. Linder 
(2004) recently pointed out that if lens models can be made 
accurate at the ~1% level then strong lensing can probe the 
dark energy equation of state in a manner complementary to 
the other probes now being discussed (such as supernovae, the 
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microwave background, and weak lensing). We have argued 
that as long as lens environments are poorly understood there 
is no way for lensing to achieve this kind of accuracy. But if 
environments receive the kind of detailed study they deserve, 
the prospects are good for lensing to participate in the new era 
of precision cosmology. 
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